spectroscopy of the stellar wind in the Cygnus X-1 system 
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Abstract 

The X-ray luminosity of black holes is produced through the accretion of material from their companion stars. Depending on 
the mass of the donor star, accretion of the material falling onto the black hole through the inner Lagrange point of the system 
or accretion by the strong stellar wind can occur Cygnus X-1 is a high mass X-ray binary system, where the black hole is 
powered by accretion of the stellar wind of its supergiant companion star HDE226868. As the companion is close to filling its 
Roche lobe, the wind is not symmetric, but strongly focused towards the black hole. Chandra-HETGS observations allow for 
an investigation of this focused stellar wind, which is essential to understand the physics of the accretion flow. We compare 
observations at the distinct orbital phases of 0.0, 0.2, 0.5 and 0.75. These correspond to different lines of sights towards the 
source, allowing us to probe the structure and the dynamics of the wind. 



Oh 
I 

o 



rn 
O 



X 



1 Introduction 

1.1 Stellar winds of O stars 

Stellar winds of early (O or B) type stars are driven by the radi- 
ation pressure of copious absorption lines present in the ultra- 
violet part of the spectrum on material in the steUar atmosphere 
Therefore the winds are very strong; common mass loss 
rates are ^lO^^Mo/year. Since primaries of high-mass X-ray 
binaries are O or early B stars [5|, which radiate in the UV, 
this radiation is strong enough to produce such a wind. Ac- 
cording to simulations of line-driven winds, perturbations are 
present and dense and cool inhomogeneities are created in the 
wind f^. Larger density, velocity, and temperature variations 
compress the gas further, creating "clumps". Current knowl- 
edge about stellar winds assumes two disjunct components of 
O star winds: cool dense clumps and hot tenuous gas. Sako 
et al. |fl6|] showed that observed spectra of X-ray binaries can 
only be explained as originating from an environment, where 
the cool and dense clumps are embedded in the photoionized 
gas. 



1.2 Cygnus X-1 

Cygnus X-1 is a binary system where the X-ray source is a 
black hole iSlil, and ~ 18 Mq Q, O 9.7 lab type star HDE 
226868 is its companion jjsl]. Stellar wind accretion plays a 
major role in the mass transfer process, because Cyg X-1 be- 
longs to the High-Mass X-ray Binaries (HMXB), which is in 
contrast to Low-mass X-ray Binaries (LMXB), where Roche 
lobe overflow is more important and accretion disk accretion 
occurs. There are strong tidal interactions in the system. More- 
over, the donor star fills ~90% of its Roche volume L5,i)]. 
Therefore the wind is not symmetric, but focused towards the 
black hole |7], such that density and mass loss rate are higher 
along the binary axis. The fact that such a high percentage of 
the Roche lobe is filled, however, means that we cannot exclude 
Roche lobe overflow taking place as well. 




Figure 1: 1.5 - 12keV light curve of Cyg X-1 obtained from 
1996 to 2010 by the RXTE-ASM. Over time, transitions from 
hard to soft state occur Compared to the time spent in the 
low-luminosity hard state (blue, < 40c/s for ASM countrate) 
Cyg X- 1 spent less time in the high-luminosity soft state (red, 
-100 c/s for ASM countrate). 



1.3 Hard and soft state of Cygnus X-1 



Black hole binaries show two principal types of emission called 
the hard or soft state, which differ in the shape of the X-ray 
spectrum, the timing properties and the radio emission. Cyg 
X-1 spends most of the time in the hard state source with a 
hard, exponentially cut-off powerlaw spectrum, strong short 
term variability and steady radio emission However, 
transitions between hard and soft states are observed (Fig.[TJ. 

How exactly the wind properties differ between states and 
what triggers state transitions are questions remaining to be an- 
swered. One possibility is that states correspond to different 
configurations of the accretion flow |17|, which cause differ- 
ences in energy dissipation. Another possibility is that changes 
of the wind properties themselves trigger state transitions ilOT . 
The mass transfer process in either case (HMXB or LMXB) 
provides extremely efficient energy release and produces lu- 
minosities of ~ lO''^ erg/s in general. For such a luminosity, 
which is typical for Cyg X-1, the X-ray source produces a con- 
siderable feedback on the wind by photoionization of its nearby 
environment [2J, which contributes to the complex wind struc- 
ture. 
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Figure 2: a) Polar view of the Cyg X-1 orbit and illustration 
of observation coverage with Chandra. Before 2010, January, 
there were 13 observations available, mostly covering the part 
of the orbit around phase <^=sO. Another observation (ObsID 
11044) was obtained at 0^0.5. A short observation (ObsID 
2742) was the only one at this phase before, but since it was 
obtained in TE mode during the soft state of the source, it 
strongly suffered from pile-up. Full lines (dashed lines) dis- 
play TE mode (CC mode) observations. Changes from blue 
to red colour correspond to changes from hard to soft state. 

b) Colorcoded orbital phases corresponding to lines of sight 
towards Cyg X-1 taking into account an inclination of ~ 35°. 

c) Highlighted observations at phases of <pxQ.O (ObsID 3814 
and ObsID 8525), 0.2 (ObsID 9847), 0.5 (ObsID 11044) and 
0.75 (ObsID 3815), which are analyzed and compared in this 
work. 

2 Observations and data analysis 
2.1 Observations and orbital coverage 

The spectra used for our analysis were obtained by the High 
Energy Transmission Grating Spectrometer (HETGS - HETG 
in combination with ACIS, Advanced CCD Imaging Spectrom- 
eter) on board the Chandra observatory. Chandra ACIS obser- 
vations are performed in two different modes: timed exposure 
(TE) mode and continuous clocking (CC) mode. In TE mode, 
the CCD is exposed for some time and then its data are trans- 
fered to the frame store, which is read out during the next expo- 
sure. The readout time required for the full frame store is 3.2s. 
In CC mode, the columns are read out continuously, which re- 
duces the readout time to 3ms f^Q. When the source is very 
bright, more than one photon may reach the same pixel in one 
frame time (pile-up). These photons are misinterpreted as one 
single event with higher energy. CC mode is usually used to 
avoid pile-up. 

High-resolution spectra of persistently bright sources like 
Cyg X-1 provide the unique possibility of probing the structure 
of the wind directly. However, this structure and therefore also 
the properties of the wind (density, velocity, ionization state) 
change with different lines of sight, which correspond to the 
different orbital phases. Thus, a good coverage of the binary 
orbit is desirable. 

Observations (whether in the hard or soft state) which are 
currently available cover part of the orbit around phase 0.0, be- 
tween phases 0.7 and 0.2, and around phase 0.5, with the latter 
only obtained in January, 2010 (Fig. |2^). We focus here on the 
comparison of observations obtained at four distinct phases of 
(p = 0.0, which is defined at the time of superior conjunction 
of the black hole (ObsID 3814 and ObsID 8525), 0.2 (ObsID 
9847), 0.5 (ObsID 1 1044) and 0.75 (ObsID 3815), see Fig. Et. 

'Chandra X-ray Center, The Chandra Proposers' Observatory Guide, 2009, 
|http : //cxc ■ harvard ■ edu/proposer/POG/| 




"Soft" (0.5-1.5 keV rate) / (1.5-3 keV rate) Softness Ratio 

Figure 3: All of these diagrams show a "soft softness ratio" 
on the X- and a "hard softness ratio" on the y-axis. Dipping 
produces a clear track in the color-color diagrams: Both colors 
harden towards the lower left comer, due to increased absorp- 
tion. However, during extreme dips, the soft color becomes 
softer again, which is likely due to partial coverage. 

The latter was obtained in CC mode, while all others were ob- 
tained in TE mode. This difference has no influence on our 
comparison. While the calibration of CC mode does not al- 
low for an adequate modelling of the whole continuum shape, 
local absorption lines, which are our primary interest, are not 
affected. 

The phase 0.0 coverage is extremely important, since due to 
the inclination of ~ 35° of the Cyg X-1 orbital plane, it cor- 
responds to looking through the densest part of the wind close 
to the stellar surface (Fig. |2j3). The distribution of X-ray dips 
with orbital phase peaks around phase 0.0 1 1]. The observation 
around phase 0.5 provides a great opportunity to close a gap in 
defining the general picture of the wind structure. 

While all recent Chandra observations caught Cyg X-1 in the 
hard state at < 100 c/s (of Chandra countrate), comparable to 
the observation at (p^Q il2ll . the spectrum was softer and the 
flux was more than twice as high during the observation at 
(p^O.l . The light curves at (p^O are modulated by strong and 
complex absorption dips, but dipping occurs akeady at <p^Q.7 
and has not ceased at <px0.2, though the dip events seem to be- 
come shorter with distance from 0=0. The light curve at 
is totally free of dips, yielding 30 ks of remarkably constant 
flux. 

2.2 Absorption dips 

According to general assumption, absorption dips, during 
which the soft X-ray flux decreases sharply, originate from in- 
homogeneities -"clumps"- present in the wind, where the ma- 
terial is of higher density and lower temperature [l6|. Ac- 
cording to the softness ratios in the color-color diagram (Fig. 
O, different stages of dipping can be classified. 

Figure |4] shows the spectrum from the observation at phase 
ipxO.l (ObsID 3815) split into "dip" and "non-dip" stages in 
the wavelength interval of the Si-region between 6 A and 7.5 A. 
While absorption lines of Si xiv and Si xm are already present in 
the non-dip spectrum, the dip spectra contain additional strong 
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Figure 4: "Dip" and "non-dip" spectra from the observation 
at phase (p^Q.l are shown for comparison of absorption Unes 
present in different stages of dipping. The reduction influx in 
the spectra is real and due to dips. In the "non-dip" spectrum, 
only absorption lines of Si xiv and Si xiii are present, while in 
the "dip" spectrum the whole series of Si xii-vii appears. 

absorption lines that can be identified with Ka transitions of 
lower ionized Si xii-vii. The strength of the low-ionization lines 
increases with the degree of dipping, indicating that the lat- 
ter is related to clumps of lower temperature. Moreover, the 
clumps are of higher density than their surroundings. In 2008 
our group organized a multi-satellite observational campaign 
with XMM-Newton,Chandra (ObsID 8525 and ObsID 9847), 
Suzaku, RXTE, INTEGRAL and Swift observing Cyg X-1 si- 
multaneously. Dips shortly after phase (p^O were so strong, that 
they were seen by all the instruments involved in the campaign, 
even RXTE-PCA or INTEGRAL-lSGRl (jj . _ 

The light curve from XMM-Newton is shown in Fig. |5] [ 1 1 ] . 
Where the dips occur in the light curve, the hydrogen column 
density, A^h, increases strongly. 

Orbilal phase 
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Day of 2008 April (MJD - 54556) 

Figure 5: Absorption dips and scattering of X-rays as seen with 
XMM-Newton, EPIC-pn. a) The light curve in the energy band 
0.3 - lOkeV shows absorption dips, which are identical to dips 
observed by Chandra shortly after (paQ. b) Hydrogen column 
density of the neutral absorption model increases strongly when 
dips occur, c) Relative flux normalization constant, which is 
consistent with the scattering trough seen in hard X-rays [11]. 

As shown in the third panel, however it is not only the pure 
absorption which causes the dips. Thomson scattering con- 
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Figure 6: Absorption and P Cygni profiles of Sixiv, Mgxii 
and Nex. While the observations at 0*0 and <f)xQ.75 show 
clear absorption profiles, although redshifted for (pxO.75, the 
observation at (f)xQ.5 shows emission at Vi-ad~0 and blueshifted 
absorption. 



tributes during these times, causing longer time scale varia- 
tions, also at hard X-rays. A possible explanation is in the ex- 
istence of dense and (nearly) neutral clumps (causing the sharp 
dips) embedded in ionized halos (causing the scattering) 111 ill . 



2.3 Spectroscopy 

We separate the "non-dip" and the "dip" parts of the obser- 
vations. The "non-dip" spectrum is extracted from the least 
absorbed phases at the upper right corner of the color-color di- 
agram (except for ObsID 3815) and the spectroscopic results 
here refer to these "non-dip" phases. The highly photoionized 
wind is detected at (p^O via numerous strong absorption lines at 
V'md^O (Fig.©. 

The lack of appreciable Doppler shifts can be explained by 
the wind flow being orthogonal to the line of sight. In con- 
trast, the recent observation (ObsID 1 1044) at (pxQ.S reveals/or 
the first time for Cyg X-1 clear P Cygni profiles with a strong 
emission component at a projected velocity Vrad~0, while the 
weak absorption components occur at a blueshift of ^500- 
1 000 km/s. If we observe the same plasma in both cases, this 
indicates that the real velocity must be small, i.e., we are prob- 
ing a dense, low-velocity wind close to the stellar surface. The 
fact that the absorption line profiles measured at (p^QJ5 are 
redshifted by x!200-300km/s indicates that the wind flow is 
not radial from the star, as a radial wind (i.e., directing away 
from the star) would always give a blueshifted velocity when 
projected onto the line of sight at phases = 0.25 - 0.75. 

3 Summary 

The new Chandra observation of Cygnus X-1 at orbital phase 
0.5 obtained in January 2010 allows us to compare observations 
at the four distinct orbital phases 0.0, 0.2, 0.5 and 0.75. With 
such a coverage, the full structure of the wind starts to reveal it- 
self. At phase 0.0 we look through the densest part of the wind, 
as it is focused towards the black hole. The light curve is mod- 
ulated by strong absorption dips. The flux decreases strongly 
during such dips, consistent with being caused by dense and 
cool clumps of material embedded in the more ionized wind. 
While absorption lines of Si xiv and Si xiii are already present 
in the non-dip spectrum, in the dip spectra also Ka transitions 
of lower ionized Si appear, whereas the strength of these lines 
increases with the degree of dipping. An especially interest- 
ing result is the totally flat light curve around phase 0.5. While 
dipping has started around phase 0.7, is the strongest around 
0.0 and still present at 0.2, it has vanished at 0.5. We therefore 
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proposed for the next observation between phases 0.25 and 0.4 
to investigate the transition between dipping and non-dipping 
phases. Spectroscopic analysis showed another interesting re- 
sult. In the spectrum at phase 0.5, clear P-Cygni profiles of 
Lyman a transitions were observed for the first time for Cyg X- 
1. We observe here strong emission components at a projected 
velocity Vrad~0 in contrast to pure absorption observed at phase 
0.0. Detailed modeUng of photoionization and wind structure 
is in progress. 
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